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Abstract 

We analyze the parameter space of the {fj, > 0, Aq = 0) CMSSM at large tan (3 with 
a small degree of non-universality originating from D-terms and Higgs-sfermion 
splitting inspired by SO (10) GUT models. The effects of such non-universalities on 
the sparticle spectrum and observables such as {g — 2)^, B{b ^ Xg'j), the SUSY 
threshold corrections to the bottom mass and QcDMh'^ are examined in detail and 
the consequences for the allowed parameter space of the model are investigated. 
We find that even small deviations to universality can result in large qualitative 
differences compared to the universal case; for certain values of the parameters, we 
find, even at low and mig, that radiative electroweak symmetry breaking fails 
as a consequence of either < or m^o < 0. We find particularly large departures 
from the mSugra case for the neutralino relic density, which is sensitive to significant 
changes in the position and shape of the resonance and a substantial increase in 
the Higgsino component of the LSP. However, we find that the corrections to the 
bottom mass are not sufficient to allow for Yukawa unification. 
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1 Introduction 



The simplest supersymmetric extension of the Standard Model, the MSSM, 
contains well over a hundred free parameters after supersymmetry breaking 
is taken into account. In the most general case, a large number of mixing an- 
gles and complex phases are present and this tends to result in predictions 
of flavour changing neutral current (FCNC) and CP-violating processes in 
excess of the current strict experimental bounds. Moreover, to analyze a pa- 
rameter space of such magnitude in any detail would require an enormous 
amount of computing power and the results would not be particularly illu- 
minating. To avoid these problems, simplifying theoretical constraints are 
usually imposed on the model in order to restrict the form of the soft su- 
persymmetry breaking mass and coupling matrices; i.e. they are taken to 
be proportional to the identity matrix in flavour space. Many analyses of 
recent years, for example [1,2,3,4,5] to cite but a few, have been based on 
the "Constrained" MSSM (CMSSM) or mSugra scenario in which a universal 
form for the mass and coupling matrices is assumed. In this model, gravity 
is presumed to be responsible for mediating the breaking of supersymme- 
try from a hidden sector of the theory sharing none of the Standard Model 
gauge interactions, to the visible sector [6,7,8,9,10,11,12,13,14,15,16,17,18]. 
With certain simplifying assumptions it is possible to construct models of 
supergravity that lead to this preferred form for the soft breaking parameters. 
However, supergravity theories by no means inevitably predict this universal- 
ity and, in any case, universality is readily violated below the supergravity 
scale by corrections deriving from, for example, renormalization group run- 
ning between the supergravity/ Planck scale and the GUT scale [19,20] or 
from the breaking of GUT [20,21,22,23,24,25,26,27,28,29] and/or family sym- 
metries [30,31,32,33,34,35,36,37,38]. 

In a previous paper [38] we explored the low energy constraints on an SO(IO) 
SUSY GUT model with an additional SU(3)f family symmetry [33,34,37] 
spontaneously broken in such a way that a phenomenologically acceptable 
set of Yukawa matrices can be obtained, accounting for the observed fermion 
masses and mixings including the neutrino sector. When both symmetries 
remain unbroken the universal form for the mass and coupling matrices for 
the gauginos and sfermions is ensured regardless of the supersymmetry break- 
ing mechanism whether it be gravity, gauge, or anomaly mediation. However, 
this sfermion universality is spoilt by D-terms arising from the breaking of 
the SU(3)f family symmetry. In this instance the mass squared of the third 
family of sfermions is split from the first two by around 20% with the sign 
of the splitting undetermined. We found that, in the case of decreased third 
family sfermion masses, the boundary of correct electroweak symmetry break- 
ing (EWSB), where ii vanishes, occurs at a substantially lower value of mo 
than for the universal case. Correspondingly, a new area of parameter space 
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allowed by the various constraints appears for large tan/3. A general con- 
clusion we reached was the allowed parameter space is very sensitive to the 
universality assumption at least for larger values of tan {3. Here we pursue this 
idea and consider the case where the dominant additional contributions to the 
sfermion mass matrices originate from the breaking of SO (10) and its sub- 
groups rather than from the breaking of the family symmetry. To isolate the 
effects of SO (10) breaking we shall ignore the effects of the family symmetry 
breaking. The impact of scalar soft mass non-universality, either deriving from 
SO (10) breaking or otherwise, on low energy phenomenology has been stud- 
ied before [39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59], but 
our analysis differs somewhat in perspective and should complement previous 
studies. In what follows, we will assume that SO(IO) is broken directly to the 
Standard Model or, in the case that there exists a secondary stage of breaking 
with an intermediate gauge group, for example the Pati-Salam group [60] : 

SO(IO) SU(4)ps X SU(2)l x SU(2)r SU(3)c x SU(2)l x U(1)y, 

or SU(5) [61]: 

SO(IO) — > SU(5) X U(l)x — > SU(3)c x SU(2)l x U(1)y, 

for simplicity we will assume that the secondary breaking scale is sufficiently 
close to the scale of SO (10) breaking that we can neglect the effects of RG 
running between those two scales. 

Additional £)-term contributions to the soft mass matrices arise from the 
reduction of rank of the gauge group, associated with a broken U(l) generator 
proportional to 2/ij -|- 3(5 — L)/2 where is the 3rd component of weak 
isospin pertaining to the SU(2)r subgroup of SO(IO) and {B — L) is the 
difference between baryon and lepton number [22,23,26,27,62]. The sfermion 
and Higgs soft masses squared each receive a contribution proportional to a 
quantity -D^, which can be positive or negative, that parametrizes the D-term 
contribution from the breaking of the GUT group. Besides the D-terms, there 
is no reason deriving from standard GUT scenarios why the EWSB Higgs 
soft masses should be related to the sfermion soft masses since they belong 
to different representations of the gauge group. Therefore we assume that the 
Higgs soft masses are independent of the sfermion masses in general, with m^j^ 
and ^ taking the common value m^Q at the GUT scale (since they are both 
contained in the 10 representation of SO (10)). Rather than allow mio and D'^ 
to vary arbitrarily, we take the ratio of GUT scale Higgs masses to sfermion 
masses, mio/mig, and similarly V = sign(£)^) ^J\D'^\/ mie as the independent 
variables. Therefore the soft SUSY breaking masses and couphngs are defined 



^ In our notation, as usual, m|^^ and the soft mass parameters for the Higgs 

fields that give mass to the down-type quarks and up-type quarks respectively. 
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by the following free parameters: 

mie, mio/mi6, mi/2, V, Aq, tan/3, sign(/x). 

Throughout this paper we will always take sign(/i) > in accordance with the 
present experimental deviation of {g — 2)^ from its Standard Model value and 
the value favoured by the branching ratio B{b Xs^f), tan/? = 50 because 
for much lower tan f3 no new allowed regions of parameter space were found 
that had not already been discovered in the CMSSM case, and Aq — since 
it doesn't have a particularly large effect on the allowed regions (except when 
it takes on large values) and also to keep the analysis simple so we can focus 
on the additional parameters introduced by the SO(IO) breaking effects^ . 

The structure of the paper is as follows: in Section 2 we review the various 
sources of non-universality present in a typical SO(IO) SUSY-GUT model 
focusing on the origin of the D-tcrm contributions to the soft SUSY breaking 
masses. In Section 3 we comment in detail on the effects on the sparticle 
spectrum of introducing the D-terms and splitting the Higgs from the sfermion 
soft masses. In particular this can have significant implications for the mass 
of the pseudo-scalar Higgs boson A° and for electroweak symmetry breaking 
which in turn has consequences for the gaugino masses and mixing angles, all 
this in addition to the mass splittings and shifts directly or indirectly induced 
in the sfermion spectrum. These changes will feed through to affect the results 
of calculations of the various observables we use to constrain the theory. We 
discuss the nature of these effects on tan /3-cnhanced amplitudes, e.g. {g — 2)^, 
in Section 4. We go on in Section 5 to give details of the numerical calculation 
of the sparticle spectrum, mixing angles and corrections to the couplings, and 
to describe the constraints on the parameter space such as the branching 
ratio of the inclusive b-decay B{b — > Xg^), the WMAP bound on the cold 
dark matter density JIcdm^^, and the discrepancy between the theoretical 
prediction of the anomalous magnetic moment of the muon and its measured 
value. In Section 6 we present the results of the constrained fit and discuss 
the main features. Finally, in Section 7 we summarize and conclude. 



2 SO(IO) Breaking and Scalar Masses 

As we touched on briefiy in the introduction, even assuming universality at 
the supergravity/Planck scale, Mp ~ 0(10^^ GeV), non-universality at the 
GUT scale, Mq ~ 0(10^^ GeV) can arise either from renormalization group 

^ Additional motivation for Aq = comes from the family symmetry. If the matter 
superfields Q, u'^j^, etc. are triplets of SU(3)f, as is the case in [37], then, at least in 
the limit of unbroken family symmetry, = 0- 
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running between the GUT scale and the Planck scale, from GUT threshold 
effects or from D-terms arising from the breaking of SO (10) or the SU(3)f 
family symmetry. In this analysis we can ignore the possibility of RGE effects 
since the full SO(IO) x SU(3)f group remains unbroken down to the GUT scale 
and hence sfermion universality is maintained. We also neglect GUT threshold 
effects originating from the discrepancy between Mq and the masses of the 
heavy particles, M//, that obtain their masses when the GUT group is broken. 
There are also potentially important corrections due to the splitting between 
the masses of the fermions, Mp, and bosons, Mb in the heavy superfields 
(see [20] for more details). This leads to corrections proportional to the log- 
arithms In and In^, and finite terms. These corrections depend on the 
details of the GUT scale particle spectrum and since we don't pin ourselves 
down to a particular SO(IO) model, and although they may be significant, we 
will ignore them for simplicity. Therefore, we choose only to include the D- 
term contributions to the sfermion masses that are expected to appear when 
this symmetry group is broken. 

Here we will outline the origin of the D-terms. For a detailed derivation, 
see [29,63]. Consider a U(l) subgroup of SO(IO), U(l)x, which is broken when 
the rank of the gauge group is reduced from 5 to 4. We assume that there 
exists two heavy scalars S(±) with opposite U(l)x charges, for example ±1, 
where in this case the charge X is given by 21 r + 3(5 — L)/2. The D-term 
contribution to the scalar potential corresponding to the U(l)x generator is 
given by 

where gi^ is the unified couphng and the U(l)x charges of the MSSM 

fields 0i. When the fields obtain VEVS, this leads to corrections to the 
mass terms for the MSSM scalar fields: 

8ml\4>,\' = x,gl({S^^)f-{S^.)f)\<P,\\ 

However, the magnitudes of the Z^-tcrms are not of the order of the large S{±^ 
VEVS as would be naively expected. A more detailed analysis shows that to 
leading order 

where frig^^^ are the soft SUSY breaking masses squared for the S scalars. 
Therefore the new D-term contributions are of order the SUSY breaking scale. 

Note that in the case of universal soft SUSY breaking masses these D-term 
contributions vanish. Since the S(^±) scalars are expected to belong to different 
representations of SO (10) (for example, potential 16h and 16h heavy Higgs 
representations) and indeed any family symmetry that may also exist at that 
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scale, they are not expected a-priori to have equal masses at the GUT scale. 
Even if they are equal at Mp this equality can easily be wiped out by differ- 
ences in the RG running down to Mq- Drawing an analogy with rn?^^ and m^^) 
a significant difference in Yukawa couplings could drive ?7t.|^^^ apart. Since we 
are otherwise ignorant of the size of D"^ we will take it to be a free parameter 
in our subsequent analysis, although we will take it to be relatively small. For 
SO (10), including the D-term contributions to the sfermion and Higgs masses, 
we obtain the well-known GUT scale inputs for the sfermion and Higgs soft 
SUSY breaking masses [22,23,26,27,62]: 



We use these GUT scale boundary conditions in our analysis, varying mio/mie 
from 0.75 to 1.25 and V from -0.4 to 0.4. These represent reasonably small 
perturbations about universality, but as we shall see, the effects on the allowed 
parameter space can be significant. 

In SO (10), a right-handed neutrino superfield completes the 16 representation. 
The presence of any neutrino Yukawa couplings and associated soft SUSY 
breaking parameters would, in principle, feed into the renormalization group 

equations for the remaining parameters above the scale of the right-handed 
neutrinos. However, we choose to ignore such effects in light of our ignorance 
of the details of the neutrino sector and our desire to keep the analysis simple. 



3 The Speirticle Spectrum 

The impact of the D-tcrms and the Higgs-sfermion mass splitting on the low 
scale phenomenology can be understood in terms of how the soft mass dif- 
ferences at the GUT scale, either directly or through their infiuence on the 
RGEs, affect EWSB and the sparticle and Higgs masses. Here we will dis- 
cuss aspects of the sparticle spectrum and mixings and EWSB. This section 
and the next arc meant to provide the reader with an intuitive guide to how 
the D-terms and Higgs-sfermion splitting feed through the different soft pa- 
rameters to affect the masses, mixings and observables. We will explain these 




^16 + 9wD'^ 
ml, - 2gl,D\ 



(1) 
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effects through various approximations vahd in certain circumstances before 
proceeding to discuss the more exact numerical results. Readers only inter- 
ested in the results may wish to skip to Section 5. We start by considering 
how the changes to the universal boundary conditions feed through the RGEs 
to create differences at the EWSB scale. 



3.1 Renormalization Group Evolution 

To begin with, we will implicitly assume that mio/mie = 1 to isolate the 
results of varying D'^. It turns out that the effect of the D-terms on the 
phenomenology is almost entirely due to their tree-level contribution to the 
boundary conditions at Mq since they almost completely cancel out in the 
RGEs [29]. We will not write down all of the soft mass RGEs here, as they 
are well known and can be found in [64] to 2-loop order. However, we will 
reproduce the RGE for third family approximation as an example 

of this cancellation and because it will be useful to refer to it later ^ . The mjj^ 
RGE is given by 

(2) 

where t — \n{Q/MG) and 

S = Tr{Ym'^) 

= m]j^ -m]j^+ ^ (m| - 2m?^ + - m| + m?^) . (3) 

generations 

Here, Y is the weak hypercharge generator. So at Mq, 

S = -Agl,D\ 

but the D-terms cancel out of the soft masses in the term proportional to 
\yt\^. This cancellation also occurs for the analogous terms in the RGEs for the 

other scalar soft masses at the one-loop level since, aside from the S term, the 
sfermion soft masses only enter into the RGEs in the following combinations 
(see for example [63]): 



^ The third family 1-loop approximation is sufficient to obtain a quafitative un- 
derstanding of the results which is our intention here. The numerical work, on the 
other hand, was carried out with the full two-loop RGEs. 
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Xt = 2\yt\^ (ml + mf^ + m|^) + 2\at\^ 

Xr = 2\yr\^ (ml^ + m?^ + m|J + 2|a,|2, (4) 

and the only remaining trace of the £)-terms in the RGEs hes within the terms 
proportional to S. Therefore, in the running from Mg to the EWSB scale, 
the D-terms only enter the RGEs multiplied by the (GUT normalised) weak 
hypcrcharge coupling gf and since wc are not considering large deviations from 
universality, we find this to be a very small effect. In the rest of this analysis we 
neglect such RG effects due to D-terms although of course we retain them in 
our numerical calculations. For more details, see [24,29]. Our main point here 
is that the splitting induced at the GUT scale by the D -terms does not change 
appreciably in magnitude as we evolve the RGEs down to the electroweak 
scale. 

Switching the D-terms off, we now allow mio/mie to vary. This time we need 
not worry about S since mio and rriiQ cancel out. However the splitting directly 
affects the Xj terms {i = t,b,T). If we take mio/mip > 1 the Xj are initially 
larger. This persists throughout the RG evolution and leads to an overall 
suppression of the 3rd family sfermion masses at low scales compared to the 
mSugra case since the effect of the Xi is to reduce the masses as we evolve 
from high to low scales. In the case of the Higgs soft masses, on the other 
hand, the tendency to be driven to lower values by the larger Xj is countered 
by the greater effect of the tree-level increase and so if we increase mio/rriiQ 
at Mq we increase mjj and mjj at the EWSB scale. 



3.2 The Heavy Higgs Masses 

The mass of the pseudo-scalar Higgs, deriving from the radiatively cor- 
rected Higgs potential in the tadpole formalism [65,66,67], is given by 

mlo = (m|^ - - M| - 7^e^i^(M|) 

COS ^f~^ 

-TZeUAA{m%) + — sin^ (3 + — cos^ (3. (5) 

Here, m^. = m^. — tj/fj where tj/fj are the corresponding tadpole contri- 
butions from loop diagrams, H^^ is the transverse part of the Z self-energy 
and H^^ is the self-energy. At large tan/5, cos 2/3 ~ — 1 and one can see 
that m^o depends overwhelmingly on the difference in the soft Higgs mass 



8 



parameters assuming that |m|^.| ^ M| 

Of importance is the RGE for rn\j^ — rn\j^, which is given by 

In general, 3Xt > SX^+Xt- due to the large top Yukawa coupling, causing 
to run to low values faster than and increasing the difference between the 
Higgs soft masses. Increasing mio/mie will increase Xi by a factor 2{w^q — 
"^i6)b«P ^G- Since 3|ytp > 3||/fep + li/rP, one would think at first sight 
that the additional contribution would tend to increase the mass difference as 
we run the parameters down to the EWSB scale compared to the CMSSM. 
However we have neglected an additional effect. In our parameter space the 
above inequality regarding the Yukawa couplings always holds because > 
suppresses yi, through the SUSY threshold corrections involving hg and x^i 
loops which contribute with opposite signs with the gluino loop dominating. 
The main contributions come from terms enhanced by tan j3 which arise from 
helicity-flipping mass insertions in the sparticle propagators. In this case both 
the gluino and chargino loops are proportional to The decrease in // caused 
by an increase in mio/mig means that the net size of the correction to y^, is 
reduced leaving us with a larger y^. It turns out that this more than counters 
the increase in Xt and for this reason we can expect the soft Higgs mass 
squared difference and thus m^o to be lighter in the case of increased mio/mie. 
As an aside we should note here that this is a small effect and, unfortunately 
from the point of view of SO(IO), we do not find Yukawa unification to better 
than about 20% in any part of the parameter space probed in this paper 
because yb is always too small for any reasonable range of mio/mie that yields 
allowed points in the parameter space. 

To illustrate the above result, Fig. 1 shows the running of the Higgs soft masses 
from Mq to Alz for the point in the parameter space mi/2 = """iie = 500 
GeV, for different values of mio/mig. We shall use this "typical" point in 
the (mi/2, mi6)-plane to compare the results of varying mio/mie and T>. It 
satisfies most of, if not all, the experimental constraints (depending on the 
exact values of T) and mio/mig) and avoids the particularly sensitive regions 
close to where EWSB fails. Also, all the approximations used in the analysis 
in this section are more or less valid in this region. Note that there is nothing 
"special" about setting mie = ?7t.i/2 since all the soft masses are renormalized 
so differently. The precise values of the sparticle and Higgs masses can be 
gleaned from Tabic 1, points A, B and C, and the values of the soft mass 
parameters are approximately given by the square of the light family sfermion 
masses. The sparticle and Higgs masses are also displayed diagrammatically in 
Fig. 1. These numerical results were produced using the method of Section 5. 
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Fig. 1. This plot shows the running of m?^^ and m'jj^ and a selection of sparti- 
de/Higgs masses for mi/2 = miQ = 500 GeV, Aq = 0, sign(fj, > 0), tan/? = 50 
and V = for (a) mio/niiQ = 1 (the CMSSM case), (b) mio/mig = 0.75 and (c) 
'^lo/^i-ie = 1.25. They correspond to points A, B and C respectively in Table 1. 

The D-terms have a much larger effect on the heavy Higgs masses than 
'^lol'^iQ since they spht the Higgs soft masses at tree level rather than through 
subtleties involving one-loop corrections. At M^, the splitting is 



This has important implications for the breaking of electroweak symmetry; 
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for large negative values of D^, m\o can be forced negative indicating that 
the electroweak symmetry has not been broken correctly; in other words, a 
solution to the potential minimisation conditions consistent with the measured 
value of Mz cannot be obtained without also having m^o < 0. Of course, 
negative does not automatically mean that this will occur since rrijj^ is 
renormalized differently from in a way that increases the difference m|^^ — 
rnjj^ . However, it does mean that m^o can be significantly smaller in this case 
than in the CMSSM. This is interesting because the region of parameter space 
where 2m^o ~ mAo, permitting resonant annihilation of LSP neutralinos via 
an S-channel A^, can occur in different parts of the parameter space as we 
shall see. Fig. 2 shows the running of the soft Higgs masses and part of the 
spectrum when the Z)-terms are non-zero. Again, the full spectrum can be 
found in Table 1 points D and E. 

Finally, in order to show the full contrast of all the effects of combining the 
D-terms with mio/mig, we show the plots in Fig. 3 (corresponding to Table 1 
points F and G). 
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Fig. 3. Similar to Fig. 1(a), hut in (a) T) = —0.4 and mio/rriiQ = 1.25, while in (b) 
T> = 0.4 and mio/jTiie = 0.75. They correspond to points F and G in Table 1. 



3.3 Electroweak Symmetry Breaking and 



The value of /i is Highly dependent on the Higgs soft masses. The minimisation 
of the Higgs potential gives the following equation [67]: 



^"h — "^ffo tan^ 3 1^1 
tan2/5-l 2^^^' 2' 



/^l = 2o . - - o^eH^z (7) 



in the same notation as previously. In the large tan f3 limit, the term involving 
m|^^ becomes irrelevant and m^^^ must be of 0{M'^) or negative at the elec- 
troweak breaking scale for an acceptable, i.e. > 0, solution to be found. 
In general, large radiative corrections due to the top Yukawa coupling achieve 
this. However, in the high mo, low mi/2 region of the CMSSM, sometimes 
referred to as the "focus-point" region [68,69,70,71] , the potency of these 



^ "Focus-point" refers to the fact that the RGEs exhibit an approximate focusing 
effect on m^j^ near the weak scale as mo is varied, keeping other parameters fixed. 
In other words the value of m^j^ is insensitive to mo and tends to be naturally of 
0{M'^). Choosing mi/2 to be very small and mo large, the soft masses and fi are 
naturally of order the EW scale in Eq. 7 in order to produce the correct result for 
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radiative corrections is reduced. This is because at small mi/2 the relative size 
of the terms chiefly responsible for driving m'jj^ negative, i.e. those involving 
multiplied by the soft squark masses and in the RGE for n'l'j^^, 
are smaller. One can understand this by realizing that a small value for mi/2 
means a small value for the gaugino mass M3 and it is terms proportional to 
g^llMsp that are the main factors involved in making the squark masses (and 
therefore Xf) large as they are RG evolved to low scales. Moreover, a large 
la^p, the square of the trilinear stop-stop-Higgs coupling, also helps to drive 
down m%^. Starting from = at Mq, at is driven negative below this scale 
by a term proportional to M3. The larger M3 is, the larger \at\ becomes and 
the stronger the effect on rrijj^ as the parameters are run towards Mz- If 
is small enough the radiative corrections may not be large enough to drive 
m|^2 sufficiently low to break the electroweak symmetry. For a given value of 
mi/2, it is generally possible to choose a large initial value of rriff^, i.e. ml, for 
which ~ at low scales. The boundary along which /i vanishes marks the 
border of correct EWSB, although before this line is reached, the LEP limits 
on the chargino mass are violated since as ^ drops below M2, the lightest 
chargino becomes ~ fi GeV. As is well known [20,39,40,72,73], by admitting 
non-universality in the Higgs sector, one can find successful EWSB in regions 
excluded in the case of universal scalar masses. 

D-terms, therefore, at the level of the input scale, can either help or hinder 
the EWSB process by increasing or decreasing rrijj^ at Mq since m^^^ = m^Q — 
2g\QD'^. Thus we may expect the = boundary to move to lower values 
of mi6 for < and to higher values for > for a given mi/2. In 
practice, however, by varying the D-terms alone (with mio/mie set to 1), 
the point m\o = is often reached before ^ = 0, depending on the exact 
location in the parameter space. On the other hand, setting the D-tcrms to 
and increasing mio/mie, we push up the value of m^^ and therefore 
towards zero at the EWSB scale towards without decreasing the value of m\o 
enormously. Therefore, at much lower values of mie than in the CMSSM, we 
find the boundary /i — where EWSB fails. 

We now look at a set of rather more interesting points in parameter space that 
illustrate some of the observations of this section — points close to where the 
breakdown of the radiative EWSB mechanism occurs. We also draw attention 
to another important fact, namely that, even relatively far from where EWSB 
fails, at low mi/2, high mig, A* becomes increasingly sensitive to changes in the 
D-terms. This can be traced back to the fact that the D-term contribution 
essentially provides an additive constant to mjj^ imposed at Mq, and, in our 
parameterization, it is proportional to rriiQ. In the small mi/2, large mie region 
of parameter space, in the CMSSM case, /i is quite small and a large positive 



Mz, resulting in very small overall fine tuning even though the sfermion masses can 
be ©(several x TeV) 
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Fig. 4. Similar to Fig. 1 hut in (a) = 500 GeV, mig = 990 GeV, mio/mig = 1 
and D = 0, in (b) — 500 GeV, mig = 990 GeV, mio/rn-ig = 1.25 and 

D = 0. These represent points U and V in Table 2. Note how light the chargino and 
neutralino spectra are in (b) compared to (a). 

value for for example, will result in a relatively large correction to m^^^ 
and therefore a relatively large increase in /x. 



To begin with, in Fig. 4(a) and (b) and Table 2 points U and V we compare 
a point in the CMSSM parameter space, mi/2 = 500 GeV, mig = 990 GeV, 
where /i is still fairly sizeable and the corresponding point with V = and 
rriiQ/miQ = 1.25 where fi is very small. Note that this occurs for a relatively 
low value of mig, much lower than the "focus point" region in the CMSSM, 
and in this region, /i is much more sensitive to the change in mio/mig than for 
our points in Figs. 1-3 and Table 1. This will be important later on when we 
come to discuss observables such as {g — 2)^ and B{b — *• Xg^y). Of particular 
importance is the lightness of the charginos and neutralinos, the lightest of 
which have important consequences for the cold dark matter relic density. 

In Fig. 5(a) and (b) and Table 2 we have points W and X. W is the CMSSM 
point at mi/2 = 813 GeV and mig = 1200 GeV. X shows the same point, but 
with V = —0.4 and mio/mig = 1.25. For X, due to the large effect of the 
negative D-terms on m^o, the heavy Higgs masses are tiny. Note that, even 
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Fig. 5. Similar to Fig. 1 hut in (a) mi/2 = 813 GeV, mig = 1200 GeF, mio/mig = 1 
and D = 0, m (b) = 813 GeF, rn-ie = 1200 GeF, mio/mig = 1-25 and 

D = —0.4. These represent points W and X in Ta6/e ^. Notice how light the heavy 
Higgses are in (b). 



though mio/miQ = 1.25, fi is still far from zero, around 300 GeV. 



The final pair of points we will consider are at mi/2 = 200 GeV, mig = 1300 
GeV, firstly for zero D-terms and mio/mig = 1 and secondly for V = 0.4 and 
mio/mi6 = 1. The results are shown in Fig. 6(a) and (b) and Table 2, points 
Y and Z, demonstrating the point made earlier: that /i is more sensitive to 
changes in (and mio/mie) at small mi/2 and large mig than for regions 
where mi/2 ~ mig, and is especially so close to the no EWSB boundary. Again, 
this fact will be of importance in our discussion of tan /3-enhanced amplitudes. 



3.4 The Sfermions 



We now consider the effects of the D-terms and mio/miQ on the sfermion 
spectrum. For a graphical representation of what follows the reader is referred 
back to Figs. 1-6. Due to the fact that they almost completely cancel out of 
the RGEs, the D-term contributions to the masses squared of the sfermions 
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Fig. 6. Similar to Fig. 1 hut in (a) mi/2 = 200 GeV, mig = 1300 GeV, mio/mig = 1 
and V = 0, in (b) mi/2 = 200 GeV, mie = 1300 GeV, mio/mie = 1 and P = 0.4. 
These represent points Y anc? Z in Ta6/e ^. T/ie D -terms affect /i about as much as 
they do m^o in this region of parameter space. 



evaluated at Msusy are numerically equivalent to their contributions to the 
soft mass squared parameters at the GUT scale, Eqs. 1, to a good approxi- 
mation. As we have discussed already, increasing mio/mig will in general lead 
to a decrease in the third family sparticle masses as a result of larger values 
of the Xi in the RGEs. However, as we shall see, by far the most significant 
effect from the point of view of phenomenology is the effect this has on fi. We 
now turn to examine the sparticle masses in detail. In all of what follows, we 
will treat the first two families of squarks and sleptons as being degenerate 
between the generations, i.e. mj^ = mg^, m^^ = mc^, rrij^^ = mc^, etc., which 
is accurate to a very good approximation. 



The Squarks 

In the case of the squarks the D-term effects are simple for the two light 
generations and the stops; we can expect the right-handed down type squarks 
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to be lighter/heavier than in the CMSSM as follows: 



~ mi 




where the soft masses on the RHS are the mSugra soft masses at Mz and gio 
is the unified gauge coupling evaluated at Mq- For the up-type squarks and 
the left-handed down squarks: 



where is one of or ni^. Since for the two light families there are 
strong radiative corrections coming from the gluino (M3) term in the RGEs 
untempered by large Yukawa terms, the squark masses are always significantly 
larger than mig and so the corrections coming from the D-terms arc relatively 
small. Prom our points in Figs. 1-6 one can see that this is a barely perceptible 
effect, in fact a few percent at most (in the case of the right-handed down 
squarks) . 

The stops receive a similar correction to the up-type squarks above: 



where i — 1,2, though it is a relatively larger effect, e.g. ~ 2 — 3% instead of 
< 1% for the first two family squarks (except for du) in Table 1, on account of 
the fact that the stop soft masses are substantially suppressed by the factors 
Xi in the relevant RGEs while the D-terms are unsuppressed. For mi/2 ^ rriiQ 
these percentages are much larger since, in this regime, M3 is small and the 
sfermion soft masses are no longer strongly enhanced at the weak scale by the 
corresponding terms in the RGEs; meanwhile, the £)-terms are relatively large, 
having been chosen proportional to rriiQ. As a result, the D-tcrm correction 
to the squark masses is relatively large. The D-terms do not significantly 
affect the stop mass mixing since they cancel out of the mixing term in the 
quadratic formula derived from the diagonalization of the mass matrix; the 
mixing remains almost entirely dependent on the trilinear coupling At which 
is insensitive to the sfermion soft masses. 

Changing mio/mie, which mainly affects 11, will not have a strong tree-level 
effect on the stop masses since we are in the large tan (3 regime where the 
//-dependent terms in the off-diagonal elements of the stop mass matrix are 
heavily suppressed. However, there will be a small decrease in the overall 
masses as noted above (around ~ 1 — 2% for the points in Table 1), due to large 
values of the in the RGEs. Again, the mixing angle is almost completely 
unaffected, being largely dependent on At which remains relatively constant. 
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For the sbottoms things are considerably more comphcated. First of all, we 
consider the average sbottom mass squared as compared to the CMSSM. From 
diagonalizing the tree-level sbottom mass matrix one obtains (ignoring EWSB 
D-terms and the bottom mass except where it is enhanced by tan /?) 



Here we have used a notation where ml = m^^'^^^^^\Mz), etc. to avoid 

ol Q33 ^ ' 

notational clutter. Therefore the average sbottom mass is lowered in the case 
of positive V and raised in the event of negative V. Due to the distribution of 
the D-terms in the sbottom mass matrix, the sbottom mass eigenstate that 
is mainly will vary a lot more than the one that is mainly 6^ if the L-R 
mixing is not too large. As far as the sbottom mass splitting is concerned one 
finds 



Amf^^ ~ y (m?^ - m?^ + Agf^D^) + 4|/^|2mf tan^ (3. (9) 

First off, we cannot say in general which of the two terms under the square 
root is the larger and cannot simply approximate the square root by the bi- 
nomial expansion. Any effect on the splitting induced by the D-term depends 
on whether m~ > m? or vice- versa, and also on the magnitude of m? — m3 

Ol Or ' ^ bL On 

relative to the magnitude of the D-term contribution. We identify two con- 
trasting regimes: 



• Large mi/2 and small mig: we find to be larger than due to the 
effect of a large value of the gaugino mass (M2) term in the RGE for m? 
which is absent from the RGE for the electroweak singlet mass squared fn^^. 
The fact that the difference between the left- and right-handed masses is 
due to, and grows with, mi/2 in this regime, and that we assume that the 
D-term grows with mig, means that the D-term is therefore relatively small 
compared to the mass splitting, even for \T>\ = 0.4, and has a much reduced 
effect. Nevertheless, in this region of parameter space, a positive will 
increase the splitting between the masses and a negative D'^ will reduce the 
splitting. 

• Small mi/2 and large mie'. we find the opposite is true — is less than 
m? due to the large Xt term in the m? RGE, a term which is again absent 
from the corresponding equation for . Here, the D-term can be relatively 
large compared to the difference between the soft masses and can easily be 
the dominant part of the first term under the square root. In this case, 
at least for large P, both signs of D'^ increase the mass sphtting though 
negative D"^ will clearly have a stronger effect. 

In certain regions of parameter space D"^ is such that it exactly cancels out 
the difference between the soft masses squared, leaving the minimum possible 



18 



splitting 

Am^^^ ~ 2|/i|mfetan/3. 
In intermediate cases things depend on the exact numbers. As it turns out, 
around the point mi/2 — miQ ~ 500, ~ and in this vicinity we cannot 
make a general connection between the sign ol V and the magnitude of the 
first term in the square root. There is an additional effect that comes into play. 
As we mentioned earlier, decreasing V results in a decrease in the size of /i, 
tending to decrease the contribution of the yU-dependent term to Arrij^^ though 
this a sub-dominant effect throughout most of the parameter space. Although 
the D-terms cause changes in the weak scale value of // a stronger effect can be 
achieved by varying mio/mie. Indeed we expect significantly smaller splitting 
of the sbottom masses for large values of this parameter where fi is suppressed. 
One can see these effects in the Figs. 1-6 and the corresponding tables. 

The sbottom mixing angle is likewise difficult to predict in the general case, 
also being sensitive to the difference between the two soft sbottom masses and 
to the off-diagonal /i dependent term. 



The Sleptons 

We begin by considering the sneutrinos and sleptons of the first two families. 
The left-handed slepton masses are given approximately by 

1 - . (10) 

whereas the right-handed slepton masses are given by 

(D) I 2 \ 

where i = 1,2 and again the soft masses on the RHS arc the mSugra masses 
at the EWSB scale and Qiq is the unified gauge coupling evaluated at Ma- 
li mi/2 is very small then mi, ~ mg^^ is of roughly equal magnitude to mig 
at the EWSB scale, at least for moderate to large values of mig. In this case 
the right- and left-handed sleptons are nearly degenerate. Otherwise m^. is 
substantially larger than both mg^- and mig, radiatively driven by the wino 
mass M2. The right-handed masses do not have a term proportional to M2 in 
the RGE and are always smaller than their left-handed counterparts in the 
CMSSM although they may still be significantly bigger than mie if mi/2 is 
large. Since the slepton masses are smaller than the squark masses as a result 
of not having an M3 term in their RGEs, the relative effect of the D-terms 
is much larger and can be as much as ~ 12% for V = 0.4 for the left-handed 
sleptons and ~ 4% for the right-handed sleptons in Table 1. Consequently, the 
right-handed charged sleptons can be heavier than the left-handed charged 



rri~.. ' ~ niz ~ m? 

^7 ^Li 
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sleptons even at moderate values of mi/2. Again, these effects can be seen 
graphically in Figs. 1-6 and numerically in the tables. 

Again, the third family is more complicated, but far more clear-cut than for 
the sbottoms. The r Yukawa coupling plays an important role by balancing 
the M2 term in the RGE for and making sure that m^^ is always pushed 
lower than mig at Mz and for low values of mig it is usually found that the 
stau is the LSP or is even tachyonic. The r sneutrino mass is given by Eq. 10 
with mi, replaced by mf^. The average stau mass is given by Eq. 8 with 
the squark soft masses replaced by the slepton soft masses. Again, increasing 
the value of will lower the average mass, with the predominantly tl mass 
eigenvalue (assuming such a distinction exists) varying more than the mainly 
right-handed one. For the mass splitting, the situation is clearer: 



Am^-^^ ~ sj (m?^ - m?^ - m^D"^) + 4|/i|2m2 tan^ /? (11) 

in the same notation as before. Note the change of sign of the D-tcrm relative 
to Eq. 9. Now, however, is always bigger than and we find that 
throughout the parameter space we explore that a positive value of always 
reduces the mass splitting (although it may swap which of the mass eigenstates 
is dominantly right- or left-handed) and a negative value will always increase 
the difference. 

For the 3rd family sleptons, although the /^.-dependent term in Eq. 11 grows 
as one decreases mio/mie, it is still very small relative to the difference in the 
soft masses squared due to the r mass squared suppression (even taking into 
account its tan^/? enhancement). In fact, since is twice as dependent on 
the factor through the RGEs as Tn\^, and Xj. is substantially smaller for 
mio/^iQ = 0.75 than in the CMSSM case, the difference "^1^ decreases. 
As a result the stau mass difference decreases as mio/niiQ decreases. The 
average masses, on the other hand, grow with decreasing mio/mig, also due 
to the smaller Xr in the RGEs. Again the reader is referred to Figs. 1-6 and 
the Tables 1 and 2. 



3.5 The Charginos and Neutralinos 

D-terms have a relatively small influence on the charginos and neutralinos. The 
soft gaugino masses Mj are almost entirely unaffected, but there is some effect 
on ji which governs the masses of the Higgsino-like charginos and neutralinos. 
We note that we can make a distinction between Higgsino-like and gaugino- 
like charginos and neutralinos when we have the hierarchy /i^ ^ M| ^ M^, 
where is the boson mass. This occurs in most of the parameter space 
that we examine and the mixing between gauginos and Higgsinos is relatively 
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small. Very roughly, in the situation where this hierarchy arises, we have 



m^o ~ Ml 
m^o ~ m ~± ~ M2 

^x°3 - ^xl - ^xt - (12) 



An increase in will cause a corresponding increase in relative to M2 
via the mechanism detailed earlier in Section 3.3. As we mentioned before, 
this is a relatively small effect. Varying mio/mie, on the other hand, creates 
substantial changes in fi and, as a result, the Higgsino-like ncutralinos and 
charginos can be significantly lighter or heavier than in the mSugra case. 
This is of vital importance for the calculation of the neutralino cosmological 
relic density where, as /i approaches Mi, the lightest neutralino may have a 
significant Higgsino component, enhancing the highly efficient Higgs-exchange 
annihilation channels. It is also important for any process with chargino and/or 
neutralino loops such as {g — 2)^, B{h — > X^^) and the SUSY correction to the 
bottom mass, 5m(,, which we have already touched on briefiy. We will discuss 
these constraints in a little more detail in later sections. 



Another matter that needs to be taken into account in order to discuss low 
energy observables is the dependence of the mixing angles in the gaugino sector 
on ji and the Mj. In the chargino sector, in the basis {W^, H^), at tree level, 
the chargino mass matrix 



_ / M2 ^MvKsin/5 
~ V\/2Miycos/3 n 

is diagonalizcd by a bi-unitary transformation with diagonalization matrices 
U and V, i.e. M.'^±^ = U* Aij(±V~^ . In the case where all the parameters are 
real, which is an assumption we make here, U and V are both orthogonal and 
parameterized by angles Or and 9l: 




^diag _ / cos6'z, sin6'z,\/ M2 \/2Mvi/ sin /5 \ / cos 6'r - sin \ 
\ — sm6L cosOl J \ \/2M[y cos /3 ji ) ysin^^ cos9r J 

(14) 

To order M^r divided by powers of soft masses (assuming that the difference 
(/i^ _ M|) > which is the case in most of the parameter space we anal- 
yse) , and in the large tan (3 regime where sin /5 ~ 1 and cos (3 ~ 0, it can easily 
be shown that 
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sin 



- Mif 



- Ml 
M|M^ 

cos R ~ 1 ^ 

Prom a quick glance, we can deduce that in this regime the chargino mixing 
angles will grow as /i decreases as a result of increasing mio/mig or decreasing 
D^. When //, approaches M2, as in the region close to where ^ vanishes at 
the boundary of EWSB for large mio/mig these approximations break down 
and this simple power series expansion is no longer applicable for the mixing 
angles. 

In the neutralino sector things are much more complicated since the mass 
matrix is 4 x 4. In the basis Hi, H^), at tree-level, 

/Ml -M^cos/5sin6'ty Af^ sin sin 6* 

M2 Mz cos /3 cos 9w —Mz sin {3 cos 

—Mz cos l3 sin 6w Mz cos jS cos 6w — 

Mz sin l3 sin 6w —Mz sin jS cos 6w -/J' 



V 



where tan^^y = is the tangent of the Weinberg angle. In the approxi- 

mation where all of the soft breaking parameters are real, this matrix can be 
diagonalized by an orthogonal matrix O, such that 



In what follows, we will be interested only in the lightest neutralino, x?. In 
the limit of large tan P, with fi"^ > Ml > M| and /x^ - Mf > M| (these 
relations again hold throughout the majority of the parameter space probed) 
an approximate solution can be found for the lightest mass eigenvalue and we 
can express roughly the components of the lightest neutralino, m^o ~ Mi, in 
the basis (5°, 1^°, H^, H^) as 



-0^0 OC 



1 



fiMz sin 9w 
M\Mz sin 9w 



(16) 



Thus, by decreasing fi relative to Mi, from an increase in mio/mig or a decrease 
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in D^, we obtain an increase in the Higgsino fraction of the lightest neutrahno. 
As we shall see, this can have an enormous effect on the neutralino relic density. 

The effects on the charginos and neutralinos of varying D"^ and mio/mie can 
be seen in Figs. 1-6 and in the tables. 



4 tan /3-enhanced Amplitudes 

We now try to draw some general conclusions about how the above changes 
in the masses and mixing angles should affect important observables such as 
the muon anomalous magnetic moment and the branching ratio B{h — > Xs^). 

It is well known that the Feynman diagrams for the quantities {g—2)^ [74,75,76], 
B{b — > Xs'~i) [77,78,79,80], and the SUSY threshold correction to the bottom 
mass [81,82,83] all contain, and in this case are dominated by, tan /^-enhanced 
contributions. This comes about when the required chirality flip for these pro- 
cesses arises from mass insertions in the gaugino or sfermion line as opposed 
to a mass insertion on the external fermion line. 



4..1 The Muon Anomalous Magnetic Moment 

We will consider in detail = . Throughout our parameter space, the 
dominant diagram is a chargino-sneutrino loop. In the absence of right-handed 

sncutrinos, the mass insertions can only occur in the gaugino line, with a /xr- 
vl-H vertex at one end and a vl-^l-W vertex at the other — see Fig. 7 — 
and the diagram is proportional to 11M2 tan {3. 

Switching now to the mass eigenstate basis for the exact 1-loop calculation of 
the chargino part of the amplitude, the tan /3 enhanced diagrams come from 
two loops, one containing Xi and one containing x^, in each case connecting 
left- to right-handed muons so that the helicity flip occurs in the chargino line. 
The part of the amplitude containing the Xi loop is proportional to a factor 
cos^ijsin^i and the part containing a loop is proportional to a factor 
— sin^RCOS^i. Throughout our parameter space we flnd these factors to be 
opposite in sign and so there is some cancellation between the loop diagrams 
involving the different chargino mass eigenstates. The charginos contribute 
factors to the amplitude that can be written in the form (compare with the 
formulae in, for example, [84]): 



23 





H~r 



+ 



^2 t^L 



Fig. 7. T/ie ian (3 -enhanced mass insertion diagram for the chargino contribution to 
{g — 2)^. The various factors show the couplings at each vertex/insertion. From the 
factor y^v sin (3 we deduce the diagram is proportional to tan fi. Contributions for 
which the helicity flip occurs on the external muon line carry a relative suppression 
of tan (3. 



In this equation, C is a positive constant with respect to varying Z^-terms or 
(at least up to neghgible corrections), and involves m^, /?, the Higgs 
VEV V and the gauge coupling §2- The function F is a phase space integral 
dependent on the relevant mass-squared ratio: 



We can project out the tan /3-enhanced part of the amplitude by setting sin (3 ~ 
1 and cos (3 ~ 0. The leading-order terms in the soft SUSY-breaking mass 
parameters are given by: 




(17) 





-3 + 4a; - a;^ - 21na; 



(19) 




(20) 
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Fig. 8. The loop function F(x), Eq. 19. 



where = m ? . And so 

^u, ^22 



c 



Here, C is a another positive constant: 



C 



ml 



167r2 



1.44 X 10"^ 



(21) 



Unfortunately, from the point of view of getting an intuitive idea of how the 
D-terms and sfermion-Higgs sphtting will affect a^, it doesn't make much 
sense to series expand the function F in terms of the mass ratios since they 
are frequently close to the critical value 1. However, F{x) is a monotonically 
decreasing function of x, and is shown in Fig. 8. Now we consider what happens 
when we increase the value of V from to 0.4 for the point mig = mi/2 = 500 
GeV. Two main effects occur. The first is that value of decreases (e.g. 
by ~ 55GeV (around 10%) for point E in table 1 as compared with point 
A). The second is that fi increases. The change is less significant for fi than 
for since /x is only affected at loop order through RGE effects (e.g. by 
around 30 GeV (around 5%) for point E compared with A). M2 remains 
basically unchanged. Consequently, the factor outside the large parentheses of 
Eq. 21, M2fi/ (m| (/i^ — M|)), grows as V increases. However, this is countered 

somewhat by what happens inside the parentheses. By increasing the mass 
ratios increase and one moves along the x-axis of Fig. 8 to the right causing 
a general decrease in the magnitude of the function in Eq. 19. One must be 
careful though because fj, will increase as a result of increasing D^, resulting 



25 



in an increase in ji^/rn^^ relative to Ml/rri^j^ therefore tending to create a 
larger difference between the two F{x) terms and cancelling out the effect 
of the smaller overall magnitude. Comparing the points E and A, the overall 
decrease narrowly wins and the difference between the loop functions F[x) is 
larger in A than in E. The external factor M2/X/ (m| (/x^ — M|)), on the other 
hand, is larger in E than in A and overall the contribution to is greater 
for V = 0.4 than for T> = 0. However, when one looks at point D for which 
V = —0.4, one can see that the cancellation between the loop function terms 
dominates the larger overall magnitude of each contribution and the difference 
between the F{x) for point D is also smaller than for point A (note, though, 
that our approximation is somewhat worse for point D than for points A and 
E due to the fact that /i is much closer in magnitude to M2). Table 3 shows 
the corresponding values for our approximation as compared to the calculated 
figures. 

We note here than in different parts of parameter space things change some- 
what. For example, for mig ^ fni/2, at least for mio/mig = 1, increasing 
has a stronger effect on a than it does on m| and the above situation is 
reversed; although this time the loop function for "D = 0.4 is much larger than 
that for T> — 0, the external factor M2/i/m| (/x^ — M|) is much smaller and 
the amphtude for T> — 0.4 is smaller than that for T> — 0. Note that in this 
region of parameter space, T> — —0.4 is excluded by the EWSB requirements 
and in any case, the above approximation would break down. 

To summarize: in general, it is somewhat complicated to understand intuitively 
how the D-terms affect a^, even though it is a relatively simple calculation as 
compared with, for example, B{b — > -^^57). It all depends on what part of the 

parameter space you arc in, where on the curve in Fig. 8 you are and also on 
the relative effects of the L'-terms on u with respect to . 

Increasing mio/mie has a more straightforward effect although it is still not 
immediately obvious what will happen from our approximation Eq. 21. Since 
an increase in this variable will tend to decrease //, while leaving the sneutrino 
mass and M2 relatively unchanged, from the above discussion it is clear that 
the difference between the loop functions will become smaller whereas the 
factor M2///m| (/i^ — M|) will grow. We find that throughout the parameter 
space, the second effect is greater and that increasing mio/mie causes to 
grow. 

4.2 B{h^Xsi) 

A similar discussion on B{h — > to that above would be rather more 

involved and not any more illuminating. We will limit ourselves to making a 
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few general remarks. 



In the Standard Model, the calculated branching ratio B{b — > Xsj) accounts 
very well for the experimentally observed value. In the MSSM one must also 
include contributions from the charged Higgs diagram, which always has the 
same sign as the Standard Model loop, and the chargino diagram. Ei- 
ther these additional contributions cannot be too large, or must cancel. In 
the positive fj, regime, the chargino amplitude has the opposite sign to the 
charged Higgs, allowing this cancellation and helping to keep the amplitude 
within the experimental bounds. However, in the CMSSM often the charged 
Higgs mass is large compared to the chargino masses at low vahics of mg and 
mi/2 and the charged Higgs contribution is not sufficient to compete with the 
chargino contribution. As a consequence the calculated value for B{b —>■ Xg'j) 
undershoots the experimental result and excludes a significant region of the 
parameter space. 

The tan /3-enhanced part of the chargino contribution to B{b — > Xs'~f) is a more 
complicated expression than the one for the {g — 2)^ chargino contribution for 
two reasons. Firstly, the approach to the calculation of B{b Xgj) is different 
— the strength of the QCD interaction necessitates operator product methods 
instead of straightforward one-loop renormalized perturbation theory and one 
has to calculate the SUSY and charged Higgs contributions to two different 
Wilson coefficients. Secondly, one also has to consider the fact that there are 
two different stop mass eigenstates propagating in the loop as opposed to the 
solitary sneutrino state, both of which contribute to tan /3-enhanced diagrams. 

When comparing the D-term- and mio/mig-corrected case to the standard 
mSugra scenario, an important factor in the discussion is the comparative 
sizes of the charged Higgs and Higgsino-like chargino masses, i.e. niAo and 
respectively. It turns out that in the b Xg'y case, unlike in the case for {g — 
2),,. wc find that the chargino amplitude always grows as /x becomes smaller, 
whether due to D-terms or mio/mig. Likewise, the charged Higgs contribution 
always increases with decreasing mj/± (or equivalcntly m^o). However, the 
question is not only to do with which of m^io or n changes most. Another factor 
is the overall size of the amphtudes. Even if the charged Higgs contribution 
becomes smaller relative to the chargino contribution and is less effective in 
cancelling it, it could be that at the same time the chargino contribution 
becomes too small to ruin the Standard Model prediction. In any case, it is 
the sum of the two amplitudes compared to the loop that is important 
rather than their relative changes in magnitude. 

If, for example, we decrease I) we cause a decrease in fi, but also a decrease 
in mH±- This effect is relatively larger for niHi than for /i over much of the 
parameter space (except for mi/2 <^ r^ie) and so one may expect the Higgs 
term to cancel the excessively large chargino contribution more successfully. 
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However, as mentioned above, the relative size of the changes is not the only 
important factor. The absohite size of the amplitudes matters too. It turns out 
that in the region mi/2 — ^le — 500 GeV, although the Higgs contribution 
grows relative to the chargino contribution with decreasing D^, they both 
increase by a roughly equal amount in absolute magnitude and their sum 
remains about the same. One can see that the B{b — > Xg^) values quoted for 
points A, D and E in Table 1 are almost equal. 

For the region mi/2 <S mie, increasing T> from to 0.4 (T> — —0.4 is ruled out 
by EWSB restrictions) affects /i just about as much as mH± and by quite a 
substantial amount. Here, the overall size of the SUSY contribution decreases; 
in addition, the charged Higgs contribution, although it decreases, cancels out 
the chargino contribution more successfully. As we shall see, this results in the 
excluded region becoming smaller as increases. 

A particularly interesting case is when is large and negative, close to the 
boundary where m^o vanishes. Here, mj^± is also very small although can 
remain several hundred GeV. In this case, one finds a narrow shver of accept- 
able parameter space where the charged Higgs contribution is large enough to 
counter the chargino amplitude even though the chargino amplitude is very 
large. We also found regions where the Higgs contribution is so large that the 
upper bound on the branching fraction is violated. 



4-3 Bottom Mass Corrections: Srrib 



Sparticle loop diagrams result in a finite threshold correction to the bottom 
mass. The most important contributions come from tan /3-enhanced terms 
which alter the relation between the (running DR) bottom mass and the 
bottom Yukawa coupling. The corrected expression can be approximated by 
(see, for example [81]) 
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l + 5mf + 5rr^t) , (22) 



where 
and 
with 



S^^ ^"-'"^^^ l(nq^,nl,mi), (23) 

*".r-^^^^^/H.m|,M% (24) 

xy\n- + yz\n^ + zx In - 

In magnitude the gluino loop tends to dominate the chargino loop due to 
the largeness of agrrig. Since At is always negative in the cases we consider. 
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the chargino loop does partly cancel the gluino loop, but overall the sparticle 
loops contribute with a positive sign in the case of positive fi {I{x, y, z) is pos- 
itive). This means that, in order to account for the experimentally observed 
value of the bottom mass, the bottom Yukawa coupling t/^ must decrease. In 
the CMSSM case, this decrease in yi, ruins the prospect of bottom-tau Yukawa 
unification. With variations in the D-tcrms and/or mio/rriie the squark masses 
change somewhat (as detailed in Section 3.4), but fi changes more. The gluino 
mass rrig is almost invariant. Both of these contributions grow with /j, times 
a loop function and as /i increases the overall contribution will tend to in- 
crease. Additionally, in most of the CMSSM parameter space, there exists a 
hierarchy where the squark masses are heavier than /x. In this case the func- 
tion / (jrij^ , m?_^, yU^^ decreases as jj, increases, helping to suppress the chargino 
contribution relative to the gluino contribution and thus enhancing the cor- 
rection to rrib and resulting in an even smaller y^. The best chance of ob- 
taining bottom-tau Yukawa unification is therefore with mio/mie — 1.25 and 
V = —0.4, but even in this case, throughout the (mi/2, mi6)-plane we are still 
far from unification — off by around 25% at best. Other authors [48,52,58] 
have found that, even allowing the D-terms and mio/mig to vary much more 
than we have, it is impossible to reconcile Yukawa unification with the posi- 
tive fjL preferred by {g — 2)n and B(6 — > X57). They have found, however, that 
if one splits the Higgs masses at the GUT scale while leaving the L>-terms 
zero, such unification is possible. One group [58] find that even then there 
are problems fulfilling the neutralino relic density constraint due to the large 
masses involved in the favoured region of parameter space. In the same pa- 
per they discuss possible solutions to this problem. However, it may be due 
to the limitations of the bottom-up procedure used in their (and indeed our) 
analysis which often fails to converge close to the region of incorrect EWSB. 
Another group [48,49], using a top-down algorithm, find no such problems and 
have identified regions of this parameter space consistent with both Yukawa 
unification and an acceptable neutralino relic density. 



5 Calculation and Constraints 

We use SOFTSUSY v.1.8.7 [85], one of several publicly available codes, to cal- 
culate the sparticle spectrum and mixings. The code has been augmented to 
include our SO(10)-inspired boundary conditions and a routine for the calcula- 
tion of the SUSY contribution to the muon anomalous magnetic moment using 
the formulae in [84] . SOFTSUSY uses a bottom-up routine in which various low 
energy observables such as Mz, fermion masses and gauge couplings are input 
as constraints in addition to the GUT scale boundary conditions. An iterative 
algorithm proceeds from an initial guess to find a set of sparticle masses and 
mixings consistent with the high and low scale constraints. We use full 2-loop 
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renormalization group equations for the gauge and Yukawa couplings and the 
fi parameter. For the soft masses we use the fuU 1-loop RGBs and include 
the 2-loop contributions in the 3rd family approximation. Full details can be 
found in [85]. A comparison between SOFTSUSY and similar programs, for ex- 
ample [86,87,88] was made in [89] and one can directly compare the codes 
online at [90]. 

For the calculation of the neutralino relic density and B{b — > Xsj) we use 
micrOMEGAs v.1.3.1 [91], linked to SOFTSUSY via an interface conforming with 
the Les Houches Accord [92] standard that contains all the relevant parameters 
from SOFTSUSY necessary for the relic density calculation. For details of these 
calculations, see [91] and the papers on which they were based [93,94,95,96,97,98,99,100,101,102]. 

In our analysis we impose the following constraints: 

• Direct searches 

The following lower limits from LEP provide the strongest constraints on 
sparticle masses from direct searches [103]: 

m^± > 103GeV mg^ > 99GeV. 

We include these lower bounds in our plots. 

• Muon anomalous magnetic moment 

We include the 2a bounds on the discrepancy between experiment and Stan- 
dard Model theory assuming the latest results of the calculation based on 
e^e~ data for the hadronic contribution [104] and the most recent data 
from the BNL E821 experiment incorporating the results from negative 
muons [105]. We use the values from [106] which include the recently recal- 
culated QED correction [107] and the most recent hadronic light-by-light 
contribution [108]. Similar values were obtained by an independent calcu- 
lation [109]. However this second paper does not take the new theoretical 
results [107,108] into account. Prom [106], 

- «r = (24-5 ± 9.0) X 10"^°, 

where = '■^"^•'^ . We use the 2a bound, 

6.5 X 10-^° < Sa^ < 42.5 x 10"^°, 

as the allowed range of the SUSY contribution. Due to the inconsistency 
between these results and those obtained by using r decay data, and taking 
into account the susceptibility to change of the measurement of the e'^e~ 
cross section [104], the {g — 2)^ constraint should perhaps be viewed more 
provisionally than the others. This is unfortunate since it is one of the most 
important, being the only one that unambiguously determines the sign of 
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• Branching Ratio B{b — > Xgj) 

The most recent world average for the branching ratio is [110] 

B{b ^ X,7)exp = (3.34 ± 0.38) x 10-^ 

while the current Standard Model theory value is [111] ^ 

B{b X,7)sM = (3.70 ± 0.30) x 10"^ 

We will use this Standard Model estimate of the theoretical error in our 
calculation as representative of the error to be expected in our calculation 
which includes both Standard Model and SUSY contributions. We do this by 
combining the experimental and theoretical errors in quadrature to obtain 
the following upper and lower bounds on the branching ratio at 2a: 

2.40 X 10"^ < B{b X,7) < 4.28 x 10"^ 

• Neutralino dark matter 

The analysis of the data from WMAP gives a best fit value for the matter 
density of the universe of Qmh^ — 0.135to;oo9 and for the baryon density, 
Qbh'^ = 0.0224 ± 0.0009 [112]. This implies that the CDM density is 

^CDM^^ = 0.11261q;q}8} 

at the 2(7 level. This can be an extremely stringent bound on the MSSM 
parameter space, especially in the case of small tan/?. However, for large 
tan (3 it is less restrictive due to the presence of the Higgs resonance, and 
much less so if we allow for a source of cold dark matter other than neu- 
tralinos such as axions, or some relic density enhancement mechanism such 
as non-thermal production of neutralinos (see [113] and references therein 
for more examples). In these instances, the lower bound on O^mh^ can be 
neglected. We plot values for which 

0.0945 < QcDuh'^ < 0.1287, 

and indicate the allowed regions if we choose to discard the lower bound. 
We also plot the locus of points for which m^o = 2171^0 marking the position 
of the resonance. 

• Lightest Higgs Mass niho 

Where appropriate, we also display the contour 

ruho = IM.lGeV, 

^ This takes into account only those results that include the improved ratio 
m^^^{mh/2)/m^°^'^ as opposed to mc°^'^/m^°'*^ in the {Xs^\{sc)v-A{cb)v-A\b) matrix 
element. For details sec [97]. 
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corresponding to the LEP bound on the hghtest SM Higgs boson [103] in 
the regions of parameter space where the hghtest MSSM Higgs boson is 
Standard Model hke, i.e. sin(/? — a) is almost exactly equal to 1, where a is 
the mixing angle relating the mass eigenstates to the gauge eigenstates in 
the CP-even neutral Higgs sector. As we shall see in the following analysis, 
there are regions of parameter space where the mass of the heavy CP-cvcn 
Higgs boson ttiho becomes comparable in mass to m/jo and sin(/3 — a) is not 
close to unity. In these regions the LEP bound on the lightest Higgs does 
not apply and we omit the contour in our plots. 
• Correct EWSB / Tachyons 

The boundary on which vanishes, marking the border of correct radia- 
tive clcctroweak symmetry breaking has been plotted. In the region where 
< a global minimum of the two loop effective Higgs potential cannot 
be found. Similarly, any regions in which m\o < 0, also signalling that the 
electroweak symmetry has not been broken correctly, have been excluded. 
Regions with tachyonic sfermions are likewise omitted. 



6 The Constrained Parameter Space 

In all of the following plots we take sign(/i) > 0, necessary to obtain the 
observed sign of {g — 2)^, and we take Aq = for simplicity, choosing to 
concentrate on the effects of varying and mio/rriiQ. Note that we do not 
demand Yukawa unification as an additional constraint. For > we were 
unable to find any regions of parameter space consistent with Yukawa unifica- 
tion. Although in principle mio / mig and can take on relatively large values 
we choose to explore the effects of fairly small deviations from universality. 
As a result there remains a close connection between the (mi/2, mo)-plane of 
the CMSSM and the (mi/2, mi6)-plane in the following analysis, i.e. for the 
most part, the sparticle spectrum will be similar in each case for correspond- 
ing points in the two planes, with possible important exceptions in the Higgs 
and chargino/neutralino sectors. Therefore we can still make a meaningful 
comparison between our SO(IO) scenario and the CMSSM. Our initial scan 
of the parameter space revealed that the most important effects only appear 
for large values of tan /3 and so we set tan p — 50 throughout. For this case, 
we will show that even small deviations from the CMSSM can lead to large 
changes in the topography of the allowed regions in the (mi/2, ?7T.i6)-plane. 

Fig. 9 shows the (mi/2, mi6)-plane in the standard CMSSM case. D-terms 
are set to zero and mio — rriiQ. One can see the usual features: the region 
where mi/2 ^150 GcV is ruled out for any value of mig due to the LEP mass 
bound on the lightest chargino; a large triangular area where mig — 0.5mi/2 
ruled out because the LSP is a stau; the quarter-egg-shaped region at small 
mi/2 excluded by the lower bound on B{b — > Xg^); the LEP lower bound on 
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Fig. 9. This plot shows contours in the {7121/2,1^1^) -plane for a variety of constraints. 
In this and subsequent plots we have Aq = 0, fj, > and tan/? = 50. In this figure 
we show the CMSSM parameter space: mio/rraig = 1 and P = 0. The blue (dark 
grey) strip and to the left of the black line at low mi/2 are excluded by the LEP 
bounds on m-^± and nif^o respectively; in the grey triangular region at small rtiiQ the 
LSP is the f; the purple (medium grey) region extending out to mi/2 ~ 500 GeV is 
ruled out by B{b ^sT); the orange (light grey) and yellow (very light grey) bands 
are the (g — 2)^ la and 2a favoured regions; the narrow crimson (darkish grey) 
curve satisfies the WMAP bounds; the dark red (dark grey) labelled line shows the 
exact position of the -resonance; the light pink (very, very light grey) region is 
allowed by WMAP if there exists another source of CDM; finally, any white regions 
are ruled out by the WMAP upper hound. 

the lightest Higgs boson mass, vahd for a Standard Model-hke Higgs; the arcs 
representing the la and 2a favoured regions for the muon anomalous magnetic 
moment at small to moderate mi/2 and mie; the A-resonance around the 
region where 2m^o = m^o and rapid annihilation can occur via an S-channel 
(or, sub- dominant ly, satisfying the upper bound on ^cdm^^; finally, 
there is the co-annihilation tail along the boundary demarking m^o = m^^ . As 
usual, the region allowed by all the constraints, if we take them all seriously, 
remains a narrow strip at fairly low mi/2 and mig, mainly dictated by the 
highly stringent WMAP bounds on i^cDM^^- 

Next, in Fig. 10 we look at the same plot, this time varying mio/mig while 
keeping P = 0. Fig. 10(a) shows mio/mie = 0.75. First of all, due to the 
resulting increase in yu, the B{h — Xs7) exclusion region has shrunk, especially 
for ^ mi/2, as would be expected from our earlier analysis in Section 4.2. 
This is because the chargino amplitude no longer over-cancels the Standard 
Model amplitude. The {g — 2)^ preferred regions follow a similar pattern (see 
Section 4.1 for details). Since M2 is more or less unaltered and in this case 
dominates the mass of the lightest chargino, the region excluded by the LEP 
bound on the Xi is almost unchanged. The stau LSP triangle is a little less 
restrictive since the Xj factors in the soft mass RGEs (see Section 3) are 
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Fig. 10. Same as Fig. 9, but showing the effects of varying mio/m-ig only; V = 0. In 
(a) m-io/m-ie = 0.75 and in (b) mio/mig = 1.25. The new grey triangular region at 
small mi/2,large uiiq in (b) is ruled out by < indicating that for this region 
there is no solution to EWSB. 

smaller and the lighter right-handed stau is slightly heavier as a result. Due 
to overall heavier masses, the contours of equal m^o move to the left and 
as a result of where they overlap with the contours of equal m^o along the 
line rriAo = 2m^o, the resonance rapid annihilation funnel is at a flatter 
angle as it emerges from the stau LSP boundary. The end result is a larger 
relic density and a smaller allowed region. 

Looking now at Fig. 10(b) where mio/mig = 1.25, we see a huge change. First 
of all, a large region appears in the upper-left of the parameter space which is 
excluded by the constraint > at the EWSB scale. Since /i is very small 
in the parameter space bordering this region, the lightest chargino becomes 
Higgsino-like and a thin strip appears along which the m + > 103 GeV bound 
is violated. The b Xg'j excluded region along with the zones of preferred {g — 
2)^ are enhanced, noticeably so at higher mig towards the edge of the EWSB 
limit where n is rapidly decreasing. Most interestingly, the enhancement of 
the A^ resonance annihilation of neutralinos has opened up the whole of the 
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parameter space preferred by {g — 2)^ and allowed by the other constraints, 
at least when one ignores the lower bound on Qcdm^-^- There arc two main 
reasons for this. The first is that the resonance, marked by the contour 
m^o = 2m^o, is steeper and exhibits a characteristic kink, thus bringing it to 
smaller values of mi/2- The kink originates from a complex interplay in the 
(mi/2, mi6)-planc near the EWSB boundary, between the difference {m'jj_^ — 
m'jj,^ (dependent on a particular combination of the RGE factors Xt and 
Xf, (discussed in Section 3.1) and which controls the evolution of m^o) and 
the absolute value of 111^^ (dependent on alone and which determines /i 
or equivalently the Higgsino component and mass of the lightest neutralino) . 
With the A'^ resonance at smaller mi/2, it is more effective at reducing JIcdm^^ 
which is proportional to m^o and therefore, for most of the parameter space 
(away from the EWSB boundary), roughly proportional to mi/2- The second 
reason is that neutrahno LSPs in the region of parameter space close to the 
boundary where |//| ~ have large Higgsino components (as can be inferred 
from Eq. 16). Since this region is at much lower mie and much closer to the A^ 
resonance than in the CMSSM case the effects are very much greater, namely 
that the neutralino LSP's coupling to the A^ and also its coupling to gauge 
bosons is greatly enhanced. Next to the m-± = 103 GeV boundary there 
is significant co-annihilation with charginos. All this conspires to produce a 
greatly diminished relic density over much of the parameter space, despite the 
fact that the GUT scale boundary conditions do not deviate massively from 
the universal case. 

In order to compare with Figs. 9 and 10, Figs. 11, 12, 13 and 14 show the 

(mi/2, '"^ii6)-plane for (a) mio/mia = 1. (b) mio/mie = 0.75 and (c) mio/mie = 
1.25, this time in order of increasing V. 

Fig. 11 shows the case oi T> = —0.4. All three graphs are fairly similar al- 
though very different from the CMSSM plot. Here, several curious things have 
happened. The area in which EWSB docs not occur is due to m\o < this 
time, {g — 2)n is slightly smaller than in the V = scenario, as expected from 
Section 4.1. Also B{b Xs'~f) is substantially less restrictive at low mig and 
the excluded region drops off sharply next to the m^o = boundary where the 
charged Higgs diagram successfully cancels out the large chargino diagram. At 
larger values of rriiQ, still along the edge of m^,, = boundary, there is a nar- 
row strip excluded by the observed value for B{b Xs^), this time because 
the Higgs component of the amplitude is ovcrcompcnsating for the chargino 
component, interfering constructively with the Standard Model contribution 
and giving too large a value. For T> — —0.4, the stau LSP excluded region is 
at its largest due to the lighter right-handed stau receiving a negative con- 
tribution due to the negative D-terra. The slope of this region is greatest for 
mio/miQ = 1.25 where Xr is at its largest, decreasing the stau mass even more. 
The A-resonance in this case has a completely different shape due to the fact 
that the Higgses are very light and beyond mi/2 ^ 800 GeV and mie ^ 450 
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Fig. 11. Same as Fig. 9, hut with T> = —0.4. In (a), mio/mig = 1, in (b) 
mio/rriiQ = 0.75 and in (c) mio/rriiQ = 1.25. The grey triangular region at small 
mi/2, large niiQ in all three graphs is ruled out this time not because < 0, but 
because m^o < 0, again indicating a failure to break EWSB properly. N.B. close to 
this region the Higgs masses are all very small. 
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Fig. 12. Same as Fig. 9, hut with V = —0.2. In (a), mio/mie = 1, in (b) 
fnw/'mm = 0.75 and in (c) mio/mie = 1.25. T/ie grey triangular region at large 
mi6 in graphs (a) and (c) is ruled out because < as for the case P = 0. 
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GeV, m^o < 2m^o everywhere. The particulars of this situation produce a 
semi-circular curve in the bottom-left of the (mi/2, 'nii6)-plane and the cor- 
responding relic density curves as shown. An enhancement of the Higgsino 
component of the LSP results in a slightly different shape for Fig. 11(c). Note, 



37 




2000 

1800 
16(]0 
1400 
> 1200 
^ 1000 

600 

iOO 
200 





/ 

^ / 

, , , / , 
200 400 «00 800 1000 1200 1400 1600 1800 2000 




200 400 600 800 1000 1200 1400 1600 ISOO 2000 









pgnp — ^ — ^ 




























/ 


V 












































y 


i 
















f 














J 


/ 11- 




















/ 








L 






/ , , 









200 400 600 800 1000 1200 1400 1600 1800 2000 



Fig. 13. Same as Fig. 9, but with V = 0.2. In (a), miQ/mi^ = 1, in (b) 
fnw/fnm = 0.75 and in (c) mio/rn-ig = 1.25. The grey triangular region at large 
mi6 in (c) is ruled out because < as for the case V = 0. 



however, that in each of (a), (b) and (c), there is no boundary where /i vanishes 
(m^o always gets there first) and there are no regions where co-annihilations 
involving charginos play a significant role. 
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Fig. 14. Same as Fig. 9, but with V = 0.4. In (a), miQ/mi^ = 1, in (b) 
fnw/fnm = 0.75 and in (c) mio/rn-ig = 1.25. The grey triangular region at large 
mi6 in (c) is ruled out because < as /or t/ie case V = 0. 



Fig. 12 shows similar plots, but with V = —0.2. This time, the situation is more 
like the situation with P = in that the no EWSB region is due to < 
rather than ruAO vanishing (although m^o is very small close to the boundary 
compared to the CMSSM). For mio/mie = 1, Fig. 12(a), unlike in the CMSSM 
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case, the D-terms help to keep above zero at the top-left-hand side of 
the plot and there is a slim region where < 0. Next to it, one can see 
at mi/2 — 300 GeV, mig — 2000 GeV a tiny sliver of parameter space ruled 
out by the B{b — > Xg'y) upper bound due to the existence of a hght charged 
Higgs. The main B{b — > ^57) exclusion and the {g — 2)^ preferred regions are 
basically indistinguishable from the CMSSM. Close to the no EWSB boundary 
at moderate to large mie there is again significant neutralino annihilation 
and there exists a region of acceptable relic density due to the large LSP 
Higgsino component. As one moves away from this region, increasing mi/2, 
the proportion of Higgsino drops off and the relic density increases above the 
WMAP upper bound. Increasing mi/2 further, one eventually reaches the A- 
resonance where again the WMAP bounds are satisfied. This results in a large, 
slanted ^/-shaped band where Qcdm^^ is perfectly adjusted to account for the 
WMAP observations. Fig. 12(b) is very similar to the case with zero D-terms, 
except that it has a larger A-resonance rapid-annihilation funnel. In Fig 12(c) 
the TTiAo = 2m^o kink is more pronounced, tending towards the semi-circular 
arc of the Fig. 11 plots. 

We note that close to the boundary where radiative EWSB fails in Figs. 9, 10, 
and 12 there is a sudden dip in the {g — 2)^ preferred region and the b ^ Xg^ 
excluded region where the Higgsino mass drops rapidly with /i, becomes nearly 
degenerate with the wino, and large cancellations occur between the partial 
amplitudes corresponding to the two mass eigenstates. For the magnetic mo- 
ment of the muon, although our approximate formulae no longer apply, one 
can infer this behaviour directly from the exact result Eqs. 17 and 18 and in- 
spection of the chargino mass matrix diagonalization equation, Eq. 14. When 
11'^ ^ M| ^ M^., then 6*^ ~ ^^/j ~ tt/A and m-± ~ m-± ~ M2 — /i and 
the contributions cancel. Of course, this cancellation is not exact, and, in any 
case, we have neglected the smuon-neutralino contributions. 

Increasing V to 0.2 and then to 0.4, we obtain the plots of Figs. 13 and 14. 
Here, and referring back to Figs. 9, 10, 11 and 12, we will remark on some 
general trends: 

• There is only a small change in B{b — > ^^7) which is most apparent at 
large mie and small mi/2, due to the increased sensitivity of /i and m^o to 
the input parameters in this region. For mio/mig = 1 or 0.75, it is at a 
maximum for "D ^ 0. This happens because, for large negative D^, a large 
contribution from a propagating H"^ results in a partial cancellation of the 
chargino contribution, thus preventing its over-destruction of the Standard 
Model amplitude, while at large positive L>^, although the charged Higgs 
diagram is comparatively negligible, the chargino loop too is much smaller 
due to the increasing masses. Both of these cases lead to a decrease in the 
excluded region. On the other hand, for mio/mie = 1.25 as one increases 
D^, the B{b — > Xg^) exclusion zone continues to increase as the = 
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boundary recedes away to smaller mi/2 and larger mie and the effect of the 
charged Higgs decreases. 

• For {g — 2)n, there are three distinct situations regarding its behaviour 
when is varied, corresponding to the different values for mio/rriiQ. For 
TTiio/miQ — 1, as grows, at high mie and small mi/2, {g — 2)^^ decreases. 
This is due to the greater effect that has on //, in this region as noted at 
the end of Section 4.1. On the other hand, at points in the parameter space 
for which mi/2 ~ m^Q, such as the mi/2 = = 500 GeV point we analysed 
in detail in Section 4.1, (g — 2)^ increases with D'^. The effect of increasing 
V at mio/mi6 = 1, therefore, is to make the muon anomalous magnetic 
moment preferred region a little lower and slightly broader. For miQ / mig = 
0.75, there is no region close to the boundary of the {g — 2)^ favoured 
zone where /x is affected to such a large extent by the D-terms since we 
are significantly further away from where EWSB fails than in the previous 
case. As a result, the {g — 2)^ preferred region increases in all parts of the 
nearby parameter space as one increases D^. Finally, for mio/wiie = 1.25, 
increasing D^, causes the = boundary to be pushed back, opening 
a larger area of parameter space. This also results in a decrease in m?^. 
Consequently, the amplitude is increased over the whole region. However, 
note the dip in {g — 2)^ as becomes small, as commented on earlier. 

• For each value of mio/mig, increasing results in a smaller f LSP region. 
This is because the lighter, right-handed stau mass steadily increases relative 
to the lightest neutralino mass since m? = m^g + gfoD"^. 

• As one increases D'^, m^o increases and the A-resonance makes a more acute 
angle with the mi/2 axis, emerging from the stau LSP excluded region at 
a larger value of mi/2- Consequently, the neutralino dark matter density is 
much increased and the region allowed by WMAP is substantially reduced 
in this part of parameter space. For V = 0.4, in Fig. 14, the characteristic 
A-resonance peak had almost disappeared. 

To make a final note on these plots, the relic density curves in the plot shown 
in Fig. 13(c) resemble those of Fig. 12(a). The difference is that the < 
region is due to an increased mio/mie in Fig. 13(c) rather than a negative 
in Fig. 12(a). 

A Cautionary Note 

It should be strongly emphasized here that the exact location and form of 
the relic density curves, the location of the EWSB excluded region and also 
of the A-resonance (on both of which the relic density strongly depends) are 
not only dependent on the variation of the mSugra parameters due to D- 
terms or Higgs-sfermion splitting or any other changes one may make, but 
also on the precise value of the top mass which governs the highly important 
parameter Xt and thus the Higgs soft masses. The value used in this paper 



41 



is mf*-""^ = 174.3 GeV, the value quoted by the PDG [103], which is now at 
odds with the recently released results from the D0 collaboration [114] (and 
see also [115]), which puts the top mass at mf^^^ = 178.0 ± 4.3 GeV. Even 
this ~ 2% change would be enough to significantly change the form of our 
plots. However, even had we used this updated value for rrit, there is still the 
question of the 4.3 GeV error. Varying rrit within this would again lead to 
very different results. Acknowledging this, we also point out that the position 
and shape of the WMAP preferred region, the A- resonance, and the = 
region are also highly dependent on the exact value of tan (3. To keep things 
under control, we used a constant value of tan (3 — 50 throughout this paper. 
However, by tweaking this value and the values we set for and mio/mig to 
some small degree we would be able to account for any difference in the top 
mass to a good approximation. As such, we believe that while qualitatively 
correct, our results for the relic density, for the exact position of the boundary 
where EWSB fails and for the exact value of the ^4*^ mass, due to their strong 
dependence on tan/? and nit, and also for the values of B{b ^s^) and 
{g — 2)^ in the region of ;U ~ or m^o ~ (although they are reasonably 
stable elsewhere) should not be taken as quantitatively correct to any degree 
of accuracy. This is not a statement of the inaccuracy of the (fairly well-tested 
and state-of-the-art) computer programs we used, merely a statement about 
the sensitivity of certain observables to the input parameters especially in 
the region of parameter space close to where the radiative EWSB mechanism 
fails. For this reason we have chosen to focus on the qualitative aspects of our 
results and their origins. 



7 Summciry and Conclusions 

In this paper we have discussed the effects that small deviations from univer- 
sality, in the form of SO(10)-inspired U(l)x -D-terms and Higgs-sfermion soft 
mass splitting, have on the parameter space of the CMSSM at large tan/3, 
with II > and Ao — 0. 

In the first part of the paper we reviewed the origin of such terms before going 
on to follow, by use of various approximations, how the additional parame- 
ters "D = sign(D^) ^ I D'^ I / rriiQ and mio/miQ feed through the rcnormalization 
group equations, electroweak symmetry breaking conditions and sparticle mass 
matrices to effect changes in the sparticle mass spectrum and mixings. 

It was noted that the D-terms are closely connected to the CP-odd ^4° Higgs 
mass and that mio/mig has a large effect on the value of fi at the electroweak 
scale. At large negative values of or large values of mio/mie one can expect 
the radiative electroweak symmetry breaking mechanism to break down due 
to m\o < in the former case and < in the latter and for this to be 
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especially important for small mi/2, large mie- 

The D-terms affect the scalars at tree-lcvcl, especially the heavy Higgs A^, 
and i?^, and the right-handed down-type squarks and left-handed sleptons, 
but largely cancel out of the RGEs. On the other hand, mio/mie affects the 
Higgs scalar soft masses at tree-level and feeds through the renormalization 
group equations through large third family Yukawa couplings, noticeably af- 
fecting the third generation of squarks and sleptons, though most strongly 
inffuencing and the Higgsino-like charginos and neutralinos. We explained 
how, contrary to first impressions, increasing mio/mig also increases m^o as 
the indirect result of bottom mass threshold corrections. Wc also gave approx- 
imate results for the mixing angles in the chargino sector and the components 
of the lightest neutralino, the Higgsino components of which can be seen to 
grow with II. 

We analysed the computed spectrum away from such regions at mi/2 = fni& = 
500 GeV for V e [-0.4,0.4] and mio/mie G [0.75, 1.25] and found deviations 
in the sparticle masses compared to the universal case resulting from additive 
effects at the GUT scale coming from the L>-terms, and from loop effects 
feeding through the RGEs as a result of the Higgs-sfermion splitting. These 
changes were found to be very small for those masses predominantly dependent 
on the gluino mass M3, for example the gluino and the squarks of the first 
two families (except the right-handed down-type squarks which receive larger 
D-teim corrections), but could be sizeable for other sparticles such as the 
third family sleptons and the Higgsino-like charginos and neutralinos. It was 
noted that close to the boundary where EWSB fails, and otherwise at regions 
with small mi/2 and large mig, the effects of D-terms and mio/mig can be 
much more important. We analyzed a set of points in these more sensitive 
regions of parameter space and observed much larger deviations in the mass 
spectrum from the universal case, especially in the heavy Higgs masses, and 
the Higgsino-like charginos and neutralinos. 

In the next section wc analysed how the D-terms and mio/mig affect the 
dominant tan /3-enhanced amplitudes of quantities such as {g — 2)^, B{b — > 
Xs^) and the SUSY threshold corrections to the bottom mass. In general the 
effects were found to be relatively small in much of the parameter space, but 
subtle, and in general larger in the small mi/2, large mie region. 

For {g — 2)n we gave a detailed quantitative approximation in Section 4.1, 
assuming that /x^— M| ^ and taking the large tan /3 limit. It was observed 
that there is a trade-off between the external factor and the loop functions 
present in the approximate formula, and in regions where the approximation 
is valid it can be quite finely balanced especially in the case of varying the 
D-terms, resulting in only small deviations from the CMSSM case. We found 
that the effect of increasing D'^ away from the sensitive regions close to l/^P = 
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and 771^0 = and to where mi/2 <^ m\Q was to increase the amphtude 
due to the decreasing muon sneutrino mass despite the contrary effect of the 
increasing Higgsino mass. For the particular point mi/2 = = 500 GeV we 
gave numerical examples, verifying our approximation. In other circumstances, 
however, often the effect of increasing D"^ was the reverse, leading to a decrease 
in (^f — 2)^. Increasing mio/mie was always found to increase {g — 2)^ as a 
result of the decrease in the Higgsino mass, though a drop-off was observed 
close to where radiative EWSB fails as a result of cancellation between the 
partial amplitudes arising from the two chargino mass eigenstates. All such 
effects were, for the most part, found to be relatively small. 

In the case of B{b Xs'~f) it was discovered that for both mio/mig = 1 and 
0.75, the size of the excluded region of B{b — > Xsj) peaks for D ~ 0. This is 
because for positive D'^ the Higgsino mass increases and therefore the chargino 
contribution to the Wilson coefficients decreases and there is no longer an over- 
cancellation of the Standard Model contribution, while for negative the 
charged Higgs contribution grows rapidly and partially cancels the chargino 
amplitude again saving the successful Standard Model result over a larger 
region of parameter space. Increasing mio/rriiQ to 1.25 in general increases 
the chargino contribution resulting in a larger excluded zone. Increasing 
in this case pushes back the no EWSB boundary, decreases the charged Higgs 
contribution and results in a larger forbidden region. 

The SUSY threshold corrections to the bottom mass were found to decrease 
with increasing mio/mie and decreasing D^, a result of increasing the chargino 
contribution with respect to the dominant gluino contribution which has the 
opposite sign. This trend is what is needed for Yukawa unification favoured 
by GUT models, but the corrections were found to be far too small to achieve 
this aim. 

Increasing mio/mig and decreasing D'^ were found to increase the area ex- 
cluded by the existence of a stau LSP due to an associated decrease in the 
mass of the lighter right-handed stau. However, the change is quite small. 

The main effect of varying mio/mig and was to cause substantial alter- 
ations to the position and shape of the curves of neutralino relic density 
required to satisfy the WMAP bounds. This was found to be due to both 
large changes in the position and shape of the A-resonance region where 
rriAo — 2m^o resulting from large corrections to the Higgs and the lightest 
neutralino mass (in the region where it changes from mainly bino to mainly 
Higgsino) , and a massively increased Higgsino component of the LSP close to 
the boundary where /i — 0, strongly affecting its couplings to the A'^ Higgs 
and gauge bosons. 

A cautionary note was made regarding the sensitivity of these results to the 
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changing value of rrit and also to tan (3. 

In conclusion we note that although corrections to universality like D-terms 
and Higgs-sfermion splitting are straightforward to implement at the input 
scale, by the time they filter through to the sparticle masses and observables 
their effects can be complex and often subtle. We have attempted to clarify 
exactly what these effects arc and how they originate. The various observables 
used to constrain low energy supersymmetric models are susceptible to change 
significantly even under small corrections to the sfermion/Higgs universality 
assumed in the CMSSM, especially observables sensitive to /i and m^o such 
as the neutralino relic density, and particularly in areas of parameter space 
close to where the EWSB mechanism breaks down. As a result a much larger 
region of parameter space than in the CMSSM becomes viable. However, from 
the point of view of a successful SO(IO) unified model, there is still the persis- 
tent problem of achieving Yukawa unification in the fi > has been 
previously noted by other authors [52,48]. We hope to address this point in 
our future work. 
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1146.9 
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1147.2 
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mi 
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1111.7 


1111.4 


1112.0 


1138.1 


1084.5 


1138.5 


1084.3 




599.5 


599.5 


599.4 


649.8 


544.5 


649.9 


544.5 
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533.9 


534.0 


510.5 


556.5 


510.5 


556.4 


mi)^ 


594.5 


594.6 


594.5 


645.3 


539.0 


645.3 


539.1 


mr 


855.9 


870.3 


837.0 


844.0 


867.7 


824.8 


881.8 




1025.5 


1036.7 


1010.2 


1013.7 


1037.0 


997.9 


1048.0 




950.1 


959.0 


937.9 


951.8 


938.0 


938.6 


946.8 




1010.7 


1023.2 


992.9 


1022.2 


1008.5 


1004.4 


1020.7 
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324.3 


328.0 


357.8 


297.3 
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539.8 


603.5 


511.8 


591.2 


527.5 




531.9 


540.3 


521.0 


588.7 


468.0 


579.1 


477.7 




393.0 


395.7 


386.1 


391.0 


394.3 


380.7 


396.4 




632.7 


685.9 


559.1 


602.7 


661.8 


527.3 


713.1 


m^o 


207.8 


207.9 


206.9 


207.5 


207.5 


206.7 


207.8 


m^o 


393.1 


395.7 


386.3 


391.1 


394.4 


381.0 


396.4 


X3 


619.7 


675.6 


539.4 


587.5 


650.4 


502.1 


703.8 


m^o 


632.3 


685.4 


559.0 


602.3 


661.3 


527.3 


712.5 


mg 


1167.4 


1167.8 


1166.9 


1167.5 


1167.2 


1167.0 


1167.6 


aSUSY X 


2.18 


2.08 


2.34 


2.04 


2.34 


2.20 


2.24 


I3{b^ Xsj) X 10^ 


2.73 


2.83 


2.57 


2.74 


2.75 


2.58 


2.84 




0.151 


0.208 


0.078 


0.074 


0.632 


0.013 


0.774 



Table 1 



Here we show the mass spectra (in GeV) for the standard points, mi/2 = = 500 
GeV. The first column, A, is the CMSSM case, B-E show the effects of changing 
one of or mio/miQ, and F and G are the combination of and mio/miQ that 
contrast the most. 
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Points of Interest 


u 


V 


W 


X 


Y 


Z 


Input Parameters 




"^1/2 


500 


500 


813 


813 


200 


200 


niiQ 


990 


990 


1200 


1200 


1300 


1300 


T) 











-0.4 





0.4 


ifiio/rriiQ 


1.0 


1.25 


1.0 


1.25 


1.0 


1.0 




581.6 


103.3 


882.2 


299.2 


293.5 


562.7 


Sparticle Masses 






116.6 


116.8 


119.3 


87.2 


113.9 


114.2 




614.2 


436.3 


846.0 


88.8 


521.2 


986.7 
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436.3 


846.0 
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521.1 


986.7 




620.1 


445.0 


850.2 


127.5 


528.0 


990.2 
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1422.0 


2007.2 


1979.9 


1346.7 


1399.1 


muii 


1395.6 


1398.1 


1960.5 


1937.4 


1347.7 


1391.2 


^4 


1421.5 


1424.1 


2008.6 


1981.4 


1349.0 


1401.3 




1388.6 


1391.1 


1947.2 


2033.5 


1347.7 


1185.5 




1040.0 


1040.2 


1308.5 


1439.1 


1300.8 


1123.8 




1006.1 


1006.1 


1235.3 


1177.6 


1300.2 


1363.1 




1037.2 


1037.4 
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1431.1 


1298.5 


1121.2 
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882.3 
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Q9Q R 
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922.7 


897.9 


1166.7 


1287.2 


1141.9 


1045.1 




914.5 


894.4 


1158.9 


1284.7 


1138.5 


933.2 


TYi -t- 


396.8 


105.2 


657.8 


312.8 


152.1 


162.5 




625.0 


422.7 


938.9 


676.3 


333.9 


598.1 


Xi 


209.2 


92.6 


347.5 


292.7 


80.8 


82.6 




396.9 


118.7 


657.9 


323.1 


152.4 


162.5 




611.3 


218.9 


927.9 


365.5 


316.6 


592.1 




624.6 


422.8 


938.8 


676.3 


332.0 


596.1 




1200.8 


1198.6 


1857.2 


1855.0 


563.1 


563.0 


^SUSY ^ 


1.17 


1.40 


0.61 


0.72 


1.17 


1.11 




2.90 


2.77 


3.37 


5.09 


2.41 


2.89 




0.857 


0.009 


0.979 


0.016 


1.364 


29.414 



Table 2 

This table shows the mass spectra ( in Ge V) for three pairs of contrasting points in 
the parameter space. Points U and V compare an unremarkable point in the CMSSM 
space (U) with the same point but with ruio/miQ = 1.25 (V) demonstrating a huge 
change in ^. Points W and X compare another unremarkable point in the CMSSM 
parameter space (W) with the same point but with mio/miQ = 1.25 and V = —0.4 
(X) where the Higgs bosons are extremely light. Finally, points Y and Z compare 
two points at large miQ and small mi/2- Point Y is the CMSSM point, while point 
Z has V = 0.4. fx is roughly as sensitive to the D-terms as m^o in this region of 
parameter space. 52 



Point 

Input Parameters 



mx/2 
V 

mio/mi6 



li{Mz) 
Mo 



m? (At2 _M|) 



10^ 



5F 

5afi X 10^ (approximation) 
5afi X 10~^ (actual value) 
(^a^ X 10^ (actual value) 



A 



D 



E 



500 

500 

1.0 



500 

500 
-0.4 
1.0 



3.06 

1.25 

0.68 

0.57 
2.54 
2.35 
2.18 



3.00 

1.39 

0.83 

0.56 
2.42 
2.24 
2.04 



500 
500 
0.4 
1.0 



592.9 561.5 622.9 
601.6 651.8 546.7 
382.8 382.8 382.7 



3.30 

1.11 

0.55 

0.56 
2.68 

2.48 
2.34 



Table 3 

This table compares the approximation to the calculated value of Sa^ for points A, 
D and E in Table 1. 
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